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ABSTRACT 


This thesis is comprised of four independent computer 
programs and their related operating instructions. Each of 
these programs focuses on a particular facet of tactical 
missile design. The topics covered include guidance and 
trajectory calculations, rocket motor propulsion siZing, 
warhead design, and aerodynamic coefficient determination. 
The programs are developed from accepted mathematical proce- 
dures and are tailored to optimize operator interaction for 
educational purposes. This thesis is intended to be uti- 
lized aS a supplement to the thesis Tactical Missile Concep- 
tual Desigqgn by D.R.Redmon, Naval Postgraduate School, Sep- 
tember 1980. 
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I. INTRODUCTION 


The programs contained in this thesis were assembled 
expressly to supplement the work of Dan Redmon in his the- 
sis, Jactical Missile Conceptual Design. Iwo of the 
programs, LPATH, the trajectory model, and LAERO1, the aéero- 
dynamic coefficients program, originated in Redmon's thesis 
and were expanded/modified for use on the Naval Postgraduate 
School's new IBM 370 computer system. The other two pro- 
grams were generated for this thesis and utilize the proce- 
dures and principles outlined by Redmon. 

The specific intention of these programs is to provide 
students of tactical missile engineering and deSign with a 
method of solving complex mathematical routines rapidly and 
interactively. Each of the programs request data which are 
likely to be used as deSign parameters for the topics con- 
cerned. The programs also allow repeated operation with 
input alteration capability, allowing the user to personally 
Optimize his design. This approach was chosen to allow stu- 
dents to understand the relationships various input paranme- 


ters have with the final solutions. 





II. TRAJECTORY MODELS 


A. DESCRIPTION AND ORIGIN OF PROGRAM 

This program applies the principles of missile guidance 
laws to the terminal phase (the last 5 to 10 seconds) of a 
missile trajectory in order to determine the maximum normal 
acceleration on the missile for a given scenario. Of the 
three general guidance law categories, pursuit guidance is 
not included in the program capability. It has been found 
that pursuit guidance invariably produces a tail-chase situ- 
ation, greatly reducing an anti-air missile's effectiveness 
against Maneuvering targets of similar speed characteris- 
tics. Line-of-sight guidance and proportional navigation 
guidance are beth options of the progran. 

Figure (II-1) shows a typical encounter geometry as 
creguired for this program. The encounter is considered to 
cccur entirely within a two- dimensional plane. No differ- 
entiation 1S required or assumed concerning the orientation 
of the encounter plane. The plane may be at any angle to 
the hcrizontal as desired by the program user. The refer- 
ence direction is an arbitrary choice by the program user. 
The angles shown are positive in value, however the progran 
does not reguire positive angles. If TAL were 150 degrees, 
1t could also Le entered as -210 degrees. The IRA tern 
represents the initial range to the target. 

Tangential velocities (air speeds) of the missile and 
the target are considered by the program to be constant 
throughout the problem. Since the program concerns itself 
With only the final moments of a trajectory, this is a réa- 
sonable consideration. farget normal accelerations, when 
specified by the user, are also held constant throughout the 


prcblem for the same reason. The missile normal 
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Program Variables 

B=LOS (line of sight angle) 

At=TAL (target alpha) 

Am=MAL (missile alpha) 

Vt=TSP (target speed) 

Vm=MSP (misSile speed) 

Rt=IRT (initial range to. 
target from missiles 
launch_point 

RM=IRM (initial range to 
missile fron 
launch point) 


a 





Rt fo Doite OLesiene 
Vm yA 
VA 


Ve Angles are shown as 
positive quantities 





Reference direction 





Figure (II-1). Encounter geometry 


ine 





acceleration, however, is a function of the guidance laws 
and will vary according to the resulting flight path 
requirements. A constraint on the encounter is that the 
missile normal acceleration must be zero at the start of the 
prcblen. 

The program "flies" both the missile and the target as 
sisple points in space with no consideration given to aero- 
dynamics. The misSile will always strike the target dead 
center when it is given the proper speed advantage for the 
encounter Since there is no provision for statistical miss 
analysis. The unbounded nature of the missile normal accel- 
eration allows the missile to turn as rapidly as necessary 
to hit the target. 

This program analyzes the given encounter by time incre- 
ment calculations. As is the case with all integrations 
conducted by incremental steps, the accuracy of the results 
is a function of the increment size. The results will 
become increaSingly accurate as the time increment is made 
Smaller. However, as the time increment decreases in size, 
the length of the output becomes increasingly longer. The 
user must balance the desire for accuracy against the amount 
of time he wishes to spend on the computer terminal. 

The primary output is a tabular listing of the missile 
and target coordinates at each time increment. The coordi- 
nate frame is established within the encounter plane with 
the abscissa oriented along the reference direction. The 
problem stops once the missile has passed its closest point 
of approach to the target. Output then provides the time to 
latercept from time of problem initiation and the maxinum 
acceleration the missile was required to endure. A Versatec 
plot of overlaying successive encounters is an optional 
output. 
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This program originated as two seperate BASIC programs 
written by Redmon {Ref. 1] for use on the HP 9830 desktop 
Calculator. It was subsequently translated into FORTRAN IV 
for use on the IBM 370 computer systen. 


Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this progran 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER]. When the program is 
not waiting for data entry but is processing, type "HX" and 
press (ENTER. Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
Simultaneously. At several points in the prograa, 
the user will be directed to "CLEAR SCREEN AND INPUT Q". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
" Aw" should appear in the lower left of the screen, 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nnnna is 
your 4-digit user number, then press [ENTER]. 

4. You will now be asked for your password. Type it in 
{the characters will not appear on the screen), then press 
(ENTER). 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 


(ENTER - 


iS 





6. Type “LINK TO xxxzP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
[ENTER]. 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |JENTER|. 
8. Type “ACCESS 192 B" and press | ENTER]. 


9. Press |ALT] and Simultaneously to clear 
screen. 

10. Type “LPATH" and press [ENTER] . 

11. Input the following data as it is requested on the 
screen by typing the data and then pressing [ENTER] . Ensure 
that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis below are the program 
variable names. 


Parameter Units Value range 

Trajectory option (TITLE) none O=Line-of-sight 
1=Proportional 

navigation 

2=Both 

Time increment (DEL) seconds Larger than .0001 
times the problen 
time 

Navigation constant (NAV) none Zao £0. b. 5 

LOS Angle (LOS) degrees 0.0 to 360.0 

Initial target ee meters Larger than the 

from the launch site (IRT) misSile range 

Target speed (TSP) n/sec Larger than 0.0 

Target normal n/sec/sec Positive is to 

acceleration (TAC) target ss Lert 

Poi tial Eee he to degrees 0.0 to 360.0 the 

line of sight (TAL) 

Missile speed (MSP) n/s ec Larger than 0.0 

Initial missile Bea ale meters Such that impact 

f£rcm launch site (IRM) occurs in less 


than 10 seconds 

12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 
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13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 


press [ENTER]. 
Switch. 


15. Turn the terminal off with the red 


C. EXAMPLE PROBLEMS 


1. Example Il-4. Jine-of-sight Non-maneuvering Cross- 
ing Target 


Not to scale 
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Launch ¥ 
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Figure (II-2). Non-maneuvering crossing target 


Rt=10000 meters 
Rm=9000 meters 
Go=30.0 degrees 
OAto=-90.0 degrees 
at=0 

Vt=200 meters/second 
Vm=800 meters/second 
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Table (II-1) is the corresponding computer output for 
this encounter. AS indicated, the missile maximum normal 
acceleration is 

am= -32.00 m/sec/sec or -3.26 g's. 
Figure (II-3) is the Versatec plot of the engagement. 

2. Exagote TI-B. P ortion Navigation Maneuverin 

Crossing Target 
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Pigure (II-4%). Maneuvering crossing target 


Rt=10000 meters 

kRm=9000 meters 

Go=30.0 degrees 

Xto=-90.0 degrees 

at=156.8 m/sec/sec (16 g's) 
Vt=200 meters/second 
Vm=800 meters/second 


Navigation constant=3.06 
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Table (II-2) is the computer output for this encounter. 
The missile maximum normal acceleration is 
am= 283.94 m/sec/sec or 28.95 g's. 
Figure (II-5) is the Versatec plot of the engagement. 


oy 


sazachi 


te 






4 tyne ehg@ene aee 
ei snldacebotae Tama 

'.2E Se Geass ae 
‘+ .8 (G27 Sereanted eae 





D. PROCEDURAL FLOWC HART 


— Input data 
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ti you want just proportional 








navagation? 


ee 


Establish initial angular 
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Calculate new angular rela- | 
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E. PROGRAM CHANGES 
1. Language Translation. 

The two BASIC pro cams contained in Reference 1 were 

translated into standard FORTRAN IV. 
2. Pxroqram Condensation. 

The two individua programs wer2 combined to forma 
Single integrated routine which allows the user *o choose 
either of the two guidance laws or both for a given encoun- 
ter. The two original projrams were meshed such that only 
the actual guidance law algorithms are separate routines, 
all input and output rout nes area now common. 

Se input and Output Facility. 

The data input instructions were modified to maxin- 
ize user facility on the I3M 370 System 3278 terminals. 
Completion of data input is now followed by data feedback 
for user verification prior to actual program execution. 

Data output has been expanded to provid2 data deliv- 
ery to three destinations; to the user's terminal for obser- 
vation, to a printer file for hardcopy duplication 2f the 
terminal display, and to a plot data file for subsequent use 
by the plot program. The jiestinations are options chosen by 
the user for each execution of the program. Up to nine 
different problems can be printed and plotted each time the 
program is entered. 
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The program can now be rerun multiple times without 
exiting and re-entering each time. The user has a choice of 
either rerunning the same problem or initiating a new prob- 
lem completely. 

4, Plot Program. 

PATHPLOT FORTRAN was developed to produce a Versatec 
plot picture of the encounter. It will produce a single 
plot each time the program is entered and will plot up to 9 
engagements in an overlaying manner. This format was 
adopted to allcw comparisons of successive input data modi- 
fications. 

5. Data Qvercapacity Check. 

If the user initiates a problem requiring more than 
9999 iterations, the program will stop. The user will be 
notified of the error and given the opportunity to rerun the 
problen. 

6. Initial Missile Acceleration. 

Initial missile acceleration was removed as a user 
input variable and established as zero. Due to the mathe- 
matical nature (il.e., no phySical constraints) of the pro-~ 
gram, any "wrong" accelerations of the missile in the ini- 
tial state were immediately "corrected" by the algorithm. 
The model is better served by providing no initial accelera- 
tions. 

7. Theta Angles. 

BOth the target and the missile taeta angles (the 
angles between the velocity vectors and the reference direc- 
tion) were removed as user input variables. The program now 
Calculates the theta angles from other input data, reducing 
redundancy and possible contradiction of input data. 

3. JZnitial Conditions Perspective. 

Originally, the missile guidance calculations 

Started at t=0. Specifically, the anchor point for the line 


ZS 





of sight solution was the point where the missile commenced 
the problem, whether or not that was on the launch site. 
This produced a situation removed from reality where the LOS 
anchor point should be at the fire control location, usually 
at or near the launch site. A Similar, though less pro- 
nounced, condition existed for the proportional navigation 
sclution. The program was modified to provide proper posi- 
ticning of the external guidance reference. AS a result of 
the modification, an additional output is the correct lead 
angle for the missile at the start of the problem. This 
angle is based on the aSsumption that neither the target nor 


the missile have maneuvered prior to time 7t=0. 
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III. WARHEAD DESIGN 





A. DESCRIPTION AND ORIGIN 

This program develops a warhead uSing the same methods 
as presented by D Redmon (Ref. 1). However, its capability 
is somewhat greater and it applies the relationships in a 
slightly different manner. This program was written for 
this thesis in FORTRAN IV for use on the IBM 370 computer 
System 3278 terminals. It is configured expressly for oper- 
ator-computer interaction. 

It starts with an initial input of data from which a 
table of fragment initial velocities is generated and pre- 
sented to the uSer. From this, the user selects values for 
fragment mass and impact velocity and another table is gen- 
erated and displayed. From this one, the desired probabil- 
ity of hit given a detonation is selected and the final 
solution is produced. At various points during the opera- 
tion of the program, the user will have opportunities to 
alter or revise selected parameters. 

The program is limited to a cylindrical warhead with 
either a solid or hollow core. The fragments: are required 
to be square in shape and are sized by the program. Figure 
(III-~1) illustrates the shape of the warhead and the loca- 
tion of various input and output quantities. 

Initially, the target altitude is used to determine the 
atmospheric density, temperature and speed of sound. These 
values are, in turn, used to calculate required initial vel- 
ccities for the fragments. The program is preloaded with 
Various values for the fragment mass and impact velocity, 
which are used to generate the table of initial velocities. 
These velocities come from the following relationship: 


V (hit) =Vifexp (-ks) J 





f 
=s, faaca 


where V(hit) is the impact velocity, Vi is the initial vel- 
ocity, s is the kill radius, m is the fragment mass, Pa is 
the atmospheric density, A is the plan area of the fragment, 
and Cd is the drag coefficient ort the fragment. 

The ballistic limit velocity is calculated for the vari- 
ous presized fragments and is provided as a reference when 
choosing an appropriate impact velocity. The ballistic 
limit velocity is that velocity at which one nalf of the 
fragments will penetrate the target's skin and the other 
half will not. The empirical relationship, developed by A. 
E. Fuhs {Ref. 6], presents the ballistic limit velocity as a 
function of the fragment size to skin thickness ratio. His 
function dealt with steel fragments impacting an aluminun 
plat2. His results were modified to qualitatively account 
for different skin materials and fragment densities. 

Next computed is the fragment spray angle and the criti- 
Cal miss distance. The spray angle 1s a function of the 
initial velocity, the detonation velocity and the warhead 
length. The critical miss distance is defined as the range 
where the fragment spray exactly covers the entire target. 
The critical miss distance is used by the program to sepa- 
rate the two functions which determine the average number of 
hits received by the target. The program assumes the target 
1s always centered within the fragment spray and aligned 
perpendicular to the central ray of the spray. 

A selection of warheads is then produced, one for each 
of a preloaded set of Pd's (probability of a hit given a 
detonation) to provide the user with a Pd versus warhead 
weight/size trade off comparison. This sizing process is 
Dased entirely on the following relationship: 


Pd=1-exp (-a) 
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where a is the average number of hits. The variable a, as 
Sshewn by Redmon [Ref. 1], is a function of the cube of the 
warhead radius. The user then chooses a desired Pd which, 
in turn, produces the final warhead sizing. 

Values that are shown in tabular form for user selection 
and input into the program are not limited to those dis- 
played. Any value in between the displayed values or con- 
pletely outside of the value range may be used. The one 
exception is that Pd can never be selected to be larger than 
-999 and may even, if forced by the program limitations, be 
creguired to be lower. Since ultimately in this algoritha, 
Pd is a function of the warhead radius, the maxinum value 
for Pd may ke reduced in order to keep the radius within the 
OCiginal input parameters. The user is notified if this 
condition occurs. 

Useful reference information for some common explosives 
and case materials is contained in the following tables. 


Table (III~1). Characteristics of common explosives 


hb AOS AVE pongo wesle lied) t ) Vd (ft/sec) 
RDX - 0596 9300. 26837. 
HMX - 0665 10230. 29934. 
FETN - 0625 9300. 21230. 
TETRYL - 0585 8200. 251596 
COMP B ~ 0607 8800. 2522 « 
OCcTOL - 0650 9500. 25350. 


Table (III-2). Densities of common case materials 


Case vitesse Density fEB/Cu. in) 
ee 


Aluminua - 100 
Uranium -688 
Patanvam e903 
Lead ~410 
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mee USER ENSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press (ENTER] . When the program is 
not waiting for data entry but is processing, type “HX" and 
press [ENTER]. Both actions will return the terminal to CHS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
simultaneously. At several points in the progran, 
the user will be directed to "CLEAR SCREEN AND INPUT OO". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
" GEA should appear in the lower left of the screen, 
press and continue. 

1. Turn the terminal on with the red Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnaP", where onnn is 
ycur 4-digit user number, then press [ENTER] . 

4. You will now be asked for your password. Type it in 
{the characters will not appear on the screen), then press 
[EytEy - 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
dees not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
(ENTER). 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
(EWrER] . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |ENTER}. 
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8. Type “ACCESS 192 B" and press [ENTER]. 


9. Press |ALT] and Simultaneously to clear 
screen. 

10. Type “LBOMB" and press [ENTER]. 

11. Input the following data as it is requested on the 
screen by typing the data and then pressing [ENTER]. Ensure 
that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis are the program variable 
names. Input the following variables as decimal values; 


Explosive density (XDEN) ee lean 
Explosive Gurney constant (GC) ft/sec 
Explosive detonation velocity (VD) Etysec 
Case material density (CDEN) PD/eus ct 
Desired kill radius (RKILL) feet 
Warhead diameter (DW) inches 


Warhead Length-to-diameter ratio (LTD) 


Target length (LT) feet 
Target mean width (WT) feet 
Target vulnerability, P(k/h) (PKH) 

Target altitude (ALT) feet 
Target skin thickness (TTH) inches 
Target skin material (MAT) 1.0=Aluminun 


2.0=Fiberglass/Kevilar 
3.0=Steel 


12. After entering the above data, you will be presented 
With an initial velocity table built around your desired 
kill radius. The initial velocities will be listed as a 
function of fragment mass and impact velocity. Also pro- 
vided will be the ballistic limit velocities for each of the 


fragment masses. Input the following parameters as decimal 


values; 
Fragment mass (IFMLB(1)) grains 
Impact velocity (VHIT(1)) ft/sec 
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13. You will now be presented with a shopping list of 
warheads developed for a range of Pd's. The warheads are 
described by the following quantities: 


Warhead weight in pounds (total weight) 
Case thickness in inches 

Core diameter in inches 

Number of fragments from the warhead 
Number of fragments hitting the target 
Edge length of the fragments in inches 


Input the following parameter as a decimal value: 
Desired probability of a hit given a detonation (PDF) 


14. If you desire to rerun the problem, or want to run a 
new problem, answer the guestions appropriately when asked 
by the terminal after the execution of the current problen. 

15. To receive the printout of your encounters, answer 
“nc” to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 

16. Upon completion of the program, type "LOGOFF" and 
press [ENTER]. 


17. Turn the terminal off with the red Sua hen. 


C. EXAMPLE PROBLEMS 

1. Example III-A 

It was desired to build a warhead which would kill a 

typical cruise missile flying at 100 feet altitude. fhe 
Warhead was reguired to have a lethal radius of 50 feet with 
a Pd of 0.98. TNT was selected for the explosive and steel 
was chosen for the case material. The diameter of the mis- 
Sile was 13.5 inches. 
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Table (III-3) outlines the input parameters. Takkle 
(IiI-4) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
Was chosen to be 130 grains and the impact velocity was 
selected to be 2500 feet per second. After the Pd table was 
displayed, 0.98 was entered as the kill probability. 

Zoeeoxkanple Iit=B 

A warhead was required which would kill a typical 
manned aircraft at 30000 feet. A lethal radius of 75 feet 
was Specified. The warhead was limited in weight to 50 
pounds and needed to have a core diameter of at least 5 
inches. The diameter of the missile was 10.0 inches. HMX 
was chosen as the explosive and depleted uranium as the case 
material. 

Table (III-5) outlines the input parameters. Table 
(III-6) shows the program output. After the initial velcoc- 
ity table was displayed at the terminal, the fragment mass 
Was chosen to ke 210 grains and the impact velocity was 
selected to be 5000 feet per second. After the Pd table was 
displayed, it was determined that the fragments were too 
large and the fragment mass was then reduced to 100 grains. 
The impact velocity was also reduced to 3000 feet per sec- 
ond. When the Pd table was redisplayed, 0.995 was chosen as 
the desired kill probability. 
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INPUT DATA FOR EXAMPLE III-A 


TABLE (IITI-3). 
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OUTPUT DATA FOR EXAMPLE III-A 


TABLE (IITI-4). 
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INPUT DATA FOR EXAMPLE III-8 
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OUTPUT DATA FOR EXAMPLE IITI-B 
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D. PROCEDURAL FLOWCHART 


Inout data | 


aa: ttt. 


Calculate characteristics of 
atmosphere at altituds 








Calculate initial velocities 
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Pees ul petcne 
velocities and fragment 
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Calculate individual warhead 
solutions for eee e!: values of 
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If pd was restricted due to 
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Determine warhead solution for 
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IV. PROPULSION MOTOR SIZING 





A. DESCRIPTION AND ORIGIN OF FROGRAM 

This program provides a method for the preliminary siz- 
ing of a solid propellant rocket motor for a boost-sustain 
trajectory of a tactical missile. The analytical method was 
developed by Redmon {Ref. 1] and was expanded with the addi- 
tion of material from Platzek [Ref. 2] and Hill [Ref. 3]. 
The program was written for this thesis in FORTRAN IV for 
use on the IBM 370 computer system. Essentially, the analy- 
sis consists of sizing the booster motor from differential 
velocity and acceleration requirements with limitations 
imposed by the physical dimensions of the missile. The 
bocster is at all times considered to be a core~burning 
motor. The sustainer motor calculations are controlled Ly 
the maximum range specified by the user and by the solution 
of the booster motor. The sustainer can be either a core- 
burning or an end~burning motor. 

The rocket motor configuration is assumed to be either 
an integral bocster-sustainer motor as shown in Figure 
(IV-1) or a staged motor as in Figure (IV-2). The booster 
and the sustainer always burn exclusively, or, in other 
words, one is not burning while the other one is. The 
nozzle half angle is specified by the user consistant with 
Space available in the missile. If a staged motor is uti- 
lized, both nozzles will have the same halr angle. 

The booster calculations start by determining the amount 
of thrust needed to boost the total weight of the missile to 
its cruise velocity at the prescribed acceleration. From 
this, the necessary amount of propellant is estimated and 
the process is iterated to account for the decreasing mass 
Situation. The chamber pressure is initially estimated by 
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Staged motor 
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Minimizing the motor weight to propellant specific impulse 
ratio, as presented by Redmon (Ref. 1]. The next step of 
the program is to size the nozzle using the now known ini- 
tial pressure ratio. Once the ideal nozzle is developed, 
the chamber pressure is raised or lowered as necessary in an 
attempt to cause the nozzle exit diameter to exactly match 
the booster diameter. However, if the integral motor option 
is used, the program will drive the nozzle as small as it 
Can without violating one of the following limits in order 
to increase the probability that the sustainer will operate 
properly. The iterative process has two limits: an abso- 
lute maximum of 2000 PSI chamber pressure [Rer. 2] anda 
Minimum of 1000 PSI if the pressure had previously been 
higher. The solution can be less than 1000 PSI if the pres- 
sure remained below that level throughout the problem. Also 
limited 1s the exit pressure to ambient pressure ratio. At 
the low end, it is limited to 0.4 to prevent flow separation 
in the nozzle. At the high end, it is limited to 2.5 to 
prevent excessive underexpansion and loss of physical real- 
ity in the program results [{Ref. 4]. The burn rate is 
initialized at a potential maximum of 1.25 inches/second 
(Ref. 2] and is allowed to decrease to arrive at a compati- 
ble burn area and web thickness combination. 

The sustainer motor, in the integral motor case, is vir- 
tually a continuation of the booster solution. The initial 
thrust requirement is determined by increasing the cruise 
speed drag to account for speed loss through maneuvers. [It 
1s also then refined for weight change if any climbing or 
diving is required to reach the target altitude. The nozzle 
is the same one as developed for the booster except that the 
throat area will be expanded as a result of the erosion 
effect. From the new area ratio, a pressure ratio is deter- 
Mined. The chamber pressure and thrust coefficient are then 
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iterated until a steady chamber pressure evolves to provide 
the reguired thrust. If at any time it drops below 125 FSI, 
the program will stop since this is considered a mininmuna 
chamber pressure for proper propellant combustion. The exit 
pressure to ambient pressure ratio remains subject to the 
Same restrictions. The burn rate starts at 0.45 inches/sec- 
ond {Ref. 2} and is decreased to provide an acceptable web 
thickness and burn area. The solution can be either an end 
burner or a core burner, depending on the burn area 
required. 

The sustainer for the staged motor is solved in essen- 
tially the same manner as the booster. The two major excep- 
tions are that it can produce an end burner if the burn area 
is small enough and the thrust required is based on the 
cruise drag and not the velocity to be gained. The chamber 
pressure is limited to an absolute maximum of 800 PSI anda 
Minimum of 250 PSI if the pressure had previously been 
higher [Ref. 2]. The same pressure ratio restrictions apply 
and the burn rate starts at 0.45 inches/second. 

Other motor-nozzle combinations can be created from the 
results of this program on a first-order approximation 
basis. Figures (IV-3) and (IV-4) illustrate two possible 
methods for separating the motor nozzles without resorting 
to staging. Although these two methods will probably have 
the same weight disadvantage that staging does, they can be 
packaged in to a Smaller volume of Space. A hint for "crea- 
tive construction" is to rerun the problem after sizing the 
motors uSing a smaller misSile diameter to force the nozzles 


to a smaller size. 
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Figure (IV-3). Separate nozzles (nonstaged motors) 
The values are obtained from the staged motor option. 





ida Onell 
= 
aaa 
<t 
ped mG 
; i 
© < 
re - 
Ay W 
—) 
Ww 
fx] = 
jaa a mA 
S |e et es 
fa pond W" 
or Ay 
© 
ae ee ~~ F&I CG 
DTS es A 
pod 
fQ 


Vv DTB-+DTS*- ~S. 


a 


ER 


NLB 


“| 


SUSTAINER 
NOZZLE 
a 
C 
n 
BOOST 
NOZZLE 


DES 










DEB-+DES4 


END VIEW 





Pigure (IV-4). Concentric nozzles (nonstaged motors) 
The values are obtained from the staged motor option. 
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Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press [ENTER] . When the program is 
not waiting for data entry but is processing, type "HX" and 
press [ENTER]. Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or “MORE...." appears in 
the status area, clear the screen by pressing and 
Simultaneously. At several points in the progran, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"GAR « should appear in the lower left of the screen, 
press and continue. 4 

1. Turn the terminal:‘on with the red |O}| switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press [ENTER]. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nonnn is 
ycur 4-digit user number, then press [ENTER]. 

G4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
[EWTER]. 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
(ENTER . 

7. YOu will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press : 
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8. Type "ACCESS 192 B" and press [ENTER]. 


9. Press [ALT] and |CLEAR| simultaneously to clear 


screen. 


10. Type "LPROP"® and press |{ENTER[. 


11. Input the following data as it 1s requested on the 


screen by typing the data and then pressing |ENTER}. 
that the data is input as either decimal or integer as spe- 


Ensure 
cified. The terms in parenthesis are the program variable 
names. 

The following variables are reguired inputs for both 
mctor option problems. 


Motor option (IMOTOR) O=integral motors 
l1=staged motors 

Launch altitude (L ALT) feet 

Launch weight (WL) pounds 

Launch velocity (V BI) feet/second 

Launch elevation angle (ELB) degrees 

Boost acceleration (A) gravities 

Cruise velocity (V BF) feet/second 

Cruise velocity drag (DRAGS) pounds 

Maximum range to target (R) hautical miles 

Maximum target altitude (TALT) feet 

Booster propellant 

speciric impulse (I SP3) seconds 

Booster propellant 

density (DENSB) pounds/cu. inch 

Bocster exhaust. 

Specific heat ratio (G B) 

Sustainer propellant 

Specific impulse (I SPS) seconds 

Sustainer propellant 

density (DENSS) pounds/cu. inch 

Sustainer exhaust. 

specific heat ratio (GS) 

Nozzle half angle (ALN) degrees 
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The following variables are required inputs for the 


integral motors option only. 


Nozzle design altitude (ALTBN) feet 


Nozzle erosion rate (E R) lunches/second 
Missile diameter (D B) inches 

Case yield strength (YIELD) PSI 

Case density (D ENSC) pounds/cu. inch 


The following variables are reguired inputs for the 


staged motors option only. 


Booster design altitude (ALTBN) feet 


Booster diameter (DB) inches 

Booster case 

yield strength (YIELD) PSI 

Booster case density (DENSC) pounds/cu. inch 
Sustainer design 

altitude (ALTSN) feet 

Sustainer diameter ( D) Laches 
Sustainer case 

yield strength (YIZLD) PSI 

Sustainer case density (DENSCS) pounds/cu. inch 


12. This program will cue the user when the input paran- 
eters have dictated a scenario which either cannot be 
achieved within reality or produce less than optinun 
reguirements on the propulsion system of the missile. They 
are not definitive and are only intended to make the user 
aware of a situation which may need correction. The follow- 
ing is a list of available cue messages with short 
definitions. 

MSUSTAINER NOT CALCULATED SINCE THE BOOSTER BURNOUT 
RANGE EXCEEDS THE DESIGN RANGE." This can result fren. 
entering an extremely short range for the missile, or it 

Can be caused by a very low acceleration reguirement. 

“SUSTAINER NOT CALCULATED; THE BOOSTER NOZZLE DESIGN 

PREVENTS SUSTAINER OPERATION. RECOMMEND STAGING OR 

INDEPENDENT NOZZLES," This occurs only when ean the 
oO 


integral motors option. Then scenario described the 
program can cause the booster nozzle to be too large to 
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maintain the proper chamber, pressures when the motor has 
shifted to suStainer operation. This usually occurs 
when a large acceleration is demanded but the thrust 
required fcr the cruise trajectory 1s small. 


"BOOSTER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEE THICKNESS." This occurs quite often and Simpl 
indicates that the burn rate was decreased from 1ts 
petential physical maximum of 1.25 inches/second. 


"THE BOOSTER NOGZLE DESIGN WAS RESTRICTED DUE TO THE 
MISSILE DIAMETER," and "THE SUSTAINER NOZZLE DESIGN WAS 
RESTRICTED DUE TO THE MISSILE DIAMETER." The nozzle was 
not able to be designed for optimum pressure ratios at 
the mid point of the boost ee oats Usually, the 
exit diameter is solved larger than the missile diameter 
and is subsequently reduced to fit. 


"BOOSTER NOZZLE DESIGN [S NOT OPTIMUM DUE TO EXCESSIVE 
BOOSTER CHAMBER PRESSURES," and "SUSTAINER NOZZLE DESIGN 
IS NOT OPTIMUM DUE TO EXCESSIVE SUSTAINER CHAMBER 
PRESSURES." If the nozzle cannot be downsized without 
peceea ng peeeeres thresholds (2000 PSI for the booster 
and 800 PSI for the sustainer), the chamber pressure 1S 
held just Eelow the pressure threshold and the nozzle 
acréa ratio wiil be adjusted to allow the nozzle to fit 
in the missile. 


"THE SUSTAINER MOTOR HAS A CORE-BURNING GRAIN." The 

reguired burn area for the Sustainer was too large tco 

Beente a end burning grain with a properly realistic 
urn rate. 


"“SUSTAINER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEB THICKNESS." This indicates the sustainer burn rate 
was lowered from a potential maximum of 0.45 inches/sec- 
ond to provide a proper web thickness. 


"THES SUSTAINER MOTOR HAS AN END BURNING GRAIN." The 
required burn area for the sustainer waS small enough to 
permit an end burning grain. The burn rate is then 
adjusted to correspond with the nonreduceable burn area. 


WREESTIMATION OF LAUNCH WEIGHT IS REQUIRED FOR THESE. 

MISSILE PERFORMANCE PARAMETERS." The scenario described 
to the program produced a motor whose weight 1s either 

merase han 75% of the total or less than 25% of the 
Otal. 


NENLARGEMENT OF DIAMETER IS RECOMMENDED DUE TO A VERY 
HIGH LENGTH-TO-DIAMETER RATIO FOR THE MOTOR." This cue 
indicates the Length-to-diameter ratio 1S greater than 
15. Other components of the missile will make the mis- 
Sile's overall length-to-diameter ratio even larger. 


13. Immediately after completion of the solution, the 


prcgram will ask if you want to receive a hardcopy printout 


of that particular solution. A "yes" answer stores that 


scluticn ina file for retrieval by the user when he fin- 


ishes with the program. 
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14. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problen. 

15. To receive the printout of your encounters, answer 
tno" to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 

16. Upon completion of the program, type "LOGOFF" and 


- Switch. 


17. Turn the terminal off with the red 


cr. EXAMPLE PROBLEMS 
1. Example IV-A. Integral motors, common nozzle 





The following parameters are input for the integral 


motor example; 


Launch altitude 35.0 feet 

Launch weight 1000.0 pounds 
Launch velocity 0.0 feet/second 
Launch angle 60.0 degrees 
Average acceleration 30.0 g's 

Cruise velocity 4000.0 feet/second 


Drag at cruise velocity 1500.0 pounds 
Maximum range 20.0 miles 
Pinal (target) altitude 50000 feet 
Booster propellant 


specific impulse 260.0 seconds 
Bocster propellant ’ 
density 0.075 lbs/cu.inch 
Booster exhaust specific 

heat ratio 1.244 
Sustainer, propellant 

specific impulse 210.0 seconds 


Sustainer propellant 


density 0.065 lbs/cu.inch 
Sustainer exhaust 

Specific heat ratio lez 0 

Nozzle half angle 20.0 degrees 
Nozzle design altitude 0.0 feet 
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OUTPUT OF EXAMPLE IV-A 


TABLE (IV-1). 
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Nozzle erosion rate 0.001 inches/second 


Missile diameter 10.0 inches 


Yield strength of 


case material 180000.0 PSI 


Density of case 
materia 


Oez2662 Lbs/cu. inch 
The solution for this problem is presented in Takkle 
(IV-1). 
2- Example IV-B. Staged motors, separate nozzles 


The following are input for the staged motor prblen: 


Launch altitude 35.0 feet 

Launch weight 2200.0 pounds 
Launch velocity 0.0 reet/second 
Launch angle 30.0 degrees 
Average acceleration 2960 99 Us 

Cruise velocity 2200.0 feet/second 
Drag at cruise velocity 1000.0 pounds 
Maximum range 50.0 miles 

Final (target) altitude 75000 feet 
pooster propellant 

specific impulse 250.0 seconds 
Bocster propellant 
density -0647 lLbs/cu.inch 
Bocster exhaust. 

specific heat ratio licezzo 

Sustainer. propellant 

specific impulse 205.0 seconds 
Sustainer propellant 
density -0625 lbs/cu.inch 


Sustainer exhaust 


specific heat ratio 1.257 

Nozzle half angle 15.0 degrees 
Booster design altitude 0.0 feet 
Booster diameter 14.5 inches 
Yield strength of 

booster case 180000.0 PSI 

Density of booster 
case material 0.2662 lbs/cu.inch 


Sustainer design altitude 
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OF EXAMPLE IV-B 


OUTPUT 


TABLE (IV-2). 
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Sustainer diameter feo ei nches 


Yield strength of 
Sustainer case 180000.0 PSI 


Density of sustainer . 
case material 0.2662 lbs/cu.inch 


The output for this example is provided in Table 


(IV-2). 
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D. PROCEDURAL FLOWCHART 


Input data 


Calculate aera de drag of 

boost phase based on cruise 

speed drag and initial and 
final speeds 


Estimate thrust reguired for 
boost and calculate booster 
burn time 





@eap aa 2a aps? «2 a= Gap «a a» a = @ apa 22 22 Gee 22 2 ee = a= a» =a 


Determine average altitude 
during boost and calculate 
ambient pressure 





Refine thrust required for 
boost and det2rmnine pec of 





| Determine nass flow rats | 
Calculate pressure ratio based 


on optimum hamber ea 
to case weight relationship 


Galcilate jehrust coefficient, 
throat diameter and exit 
diameter 


Is exit diameter smaller than 
base diameter? 


NO 


| Increase chamber pressure 


ed nc eae ee eeS 
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NO Is chamber pressure larger 
eee 2 
yes 
Lower chamber pressure? 
Increase 2x1t pressure 
Is (exit pressure/ambient 
. pressure) greater than 0.4? 
NO 
Is (exit pressure/ambient 
pressure) less than 2.5? 
NO 


Set exit diameter equal to 
base diameter 














Is the nozzl2 2xit diameter 
larger than 95% of the base 
dianeter? 









Is the chamber pressure less 
than 100 PSI? 
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Set burn rate egual to 1.25 
and calculate] resulting web 
thixkness 


Determine burn area and 
required web thickness 


Are the two web thickness 
values egual? 


Calculate di nensions of the 
booster and nozzle 








Calculate required sustainer 
urn time 
Determine average altitude 


during cruis= and calculate 
ambient pressure 





Calculate thrust required to 
offset cruise drag and 
Calculate weight of sustainer 
propellant 
Determine mass flow rate and 
nozzle area ratio 








Calculate Mach of the exhaust 


eoncenene tach goa ons Sba22et. 
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Calculate nozzle pressure 
ratio and estimated thrust 
coefficient 


Calculate sustainer chamber 
pressure 
Does thrust coefficient equ 
estimated thrust coefficien 
Set them equal 











“~ Is the Durn area less than the 
motor cross sectional ar2:3? 


Is the web thickness less tha 
86% of the motor radius? 


> Yes 


Set burn ar a equal to motor 
cross-section area and lower 
burn rate to sorrespond 


Establish dimensions of sus- 
tainer motor 











Determine average altitude 
during cruise and calculate 
ambient pressure 


Calculate thrust required to 
offset cruise drag and 
calculate weight of sustainer 
propellant 


Calculate pressure ratio based 


on optimum, hamber pressure 
to case weight relationship 


Calculate thrast coefficient, 
throat dianeter and exit 
diameter 


s exit diameter smaller than 
ase diameter? 


Is chamber pressure larger 
than PLIM? 


Lower chamber pressure 
Increase 2xit pressure 


Is (exit pressure/ambient 
pressure) greater than 0.4? 





NO 
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Is (exit poieeeee ambient 
¢ pressure) Less than 2.5? 


NO 


Set @2xit diameter equal to 
ase diameter 







Is nozzle exit diameter larger 

than 95% of the base diameter? 

Is the chamber pressure less 
Ehagezou PST? 


Is burn area smaller than 
motor cross-section area? 


Calculate required web thick- 
ness and web thickness 
resulting fron the burn rate 


= Font pelea Mi poasa sees cate. 
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@mepaenpen eee Baepweeawe ae 82 2 SS Bee 2 a a2 a = 28 ee op 2 oF oe 


Set burn area equal to motor 
cross-section area and 
calculate corres ponding burn 
rate 


Calculate dinensions of sus- 
tainer motor and nozzle = 


Zero 


Print output | 


naman SEES 


[exo 
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Wo WAP RODYNANITC CORFEILCLENTS 


A. DESCRIPTION AND ORIGIN 

This program is the current edition of a program which 
originated at the Naval Ship Research and Development Center 
in 1971 (Ref. 5]. The program was written in FORTRAN II for 
use on the IBM 7090 digital computer. It was developed as a 
method for predicting the static, longitudenal aerodynamic 
characteristics of typical missile configurations with the 
ccntrcl surfaces in a "plus" attitude. The original program 
computed the aerodynamic characteristics for missiles oper- 
ating at angles of attack up to 180 degrees. The effects of 
control surface deflections for all modes of aerodynamic 
ccntrcl are taken into account. The method was based on the 
then well Known linear, nonlinear crossflow and slender kody 
theories with empirical modificaticns to provide the high 
angle of attack capability. 

The program was modified and, presumably, updated in 
1574 by F. A. Kuster, Jr., of the Naval Air Development Cen- 
ter. In 1980, the program was modified for use on the Naval 
Postgraduate School IBM 360 computer system by D. Redmon 
{[Ref. 1]. The current version of the program was modified 
fcr use on the new Naval Postgraduate School IBM 370 compu- 
ter system. It has been expanded to provide graphical pre- 
sentation of the output data. 

It must be emphasized at this point that the current 
program edition is not in a completely finished state. 
Scmewhere in the history of the program after its initial 
establishment on the IBM 360, errors were introduced during 
the modifications. At present, these errors do not prevent 
the use of the program and the output data is considered to 


be correct for trend-observance purposes. 
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Specifically, the program does not produce any drag-rise 
phenomena for either the wings or the tails when Cd is 
observed as a function of Mach number. Additionally, the 
decline of the drag coefficient above Mach 1.0 is not smooth 
or as prolonged as is found experimentally. It is very fro- 
bable that these two failings of the program are linked to a 
ccmmon error inserted accidentally in the process of tailor- 
ing the program for use on the IBM 370. In order to 
temporarily smcoth over the graphical discontinuities, expo- 
nential decay functions were inserted. They are clearly 
marked in the the program listing for removal when the pro- 
gram is corrected. 

The input to the program is composed of a detailed list- 
ing of the dimensions of the missile to be analyzed. The 
current version of the program will consider a missile which 
has four symmetrical wings and four symmetrical tails. The 
missile may be either canard or tail equipped and either 
Wing or canard or tail controlled. The program assumes the 
ccntrol surface is the tail, however, the input data is 
"mislabelled"™ to produce the proper configuration. For ins- 
tance, if the missile is a wing control missile, the wing 
data is input as the tail and the tail data as the wing. 

For a canard controlled missile, the canard data is input as 
the tail. Figures (V¥-1) and (V-3) show two typical missile 
ccnfigurations and where the input parameters for the pro- 
gram are obtained. 


Be. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
Waiting for data entry, press [ENTER] . @hen the program is 
not waiting for data entry but is processing, type "HX" and 
press [ENTER]. Both actions will return the terminal to CMS 
mode. 
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When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing and 
simultaneously. At several points in the progran, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 
This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"GEA" should appear in the lower left of the screen, 
press and continue. ] 

mee cucn the terminal on with the red . Switch. 

2. When the large "NPS" appears after about 30 seconds, 
press [RESET], then press {ENTER}. 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnonnP", where nnnn is 
ycur 4-digit user number, then press (ENTER] . 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 
(eurea] 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH", then press 
(ENTER). 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 
[eurey}. 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 


then press |[sNTER]. 
8. Type “ACCESS 192 B" and press [ENTER]. 
ae ~PLeSS and {CLEAR| simultaneously to clear 


screen. 


10. Type "LAERO1" and press [ENTER]. 


11. Input the following data as it is requested on the 
screen by typing the data and then pressing |{ENTER]. Ensure 
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that the data is input as either decimal or integer as sfe- 
cified. The terms in parenthesis below are the program 


variable names. 


Input the following as integer values unless otherwise 
noted. The integers must be two digit integers (three=03). 


Centrol constant (ICSC) Q1-Tail control 
; Q2-Wing control 
0O3-Canard control 
Nose shape (INOSE) Q1-ELlipsoid 
g4r23z3" 


Number of control deflections (IDT) Less than 11. 
You will now be asked for the control deflections in 
degrees, decimal values. 


Number of Mach numbers (IM) Less than 25. 
You will now be asked for the Mach numbers, decimal 
values. Each Mach number will produce a separate 
table and plot of output data. 


Number of angles of attack (IAL Less than 25. 
you will now be asked for the angles of attack in 
degrees, decimal values. 


Number of configurations ({NBODY No restrictions 
Each configtration will restart the program. Only 
the last configuration will produce the written 


output. 
Wing planform (ISWPW) Q1-Not delta 
Q2-Delta 
Wing position (IAFBW) OO an ae after wing. 
Q1-No body after wing 
Wing sweep constant (ISWEPW) 00-Delta planform or 
See leading 
edge 
ese leading edge 
Number of wings (NWING) O4 
Tail planform (ISWPT) O1-Not delta 
Q02-Delta 
Tail position (IAFBT) QQ-No Boy after tail 
Q1-Body after tail 
Tail sweep constant (ISWEPT) 00-Delta planform or 


Unswept leading 
edge 
Q01-Swept leading edge 
Number of tails (NTAIL) O4 . 
Input the following values as decimal numbers: 
Wing tip-to-chord ratio (XLAMW) 
Wing leading edge sweep (CLAMW) Degrees 
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Wing span including body (BX) 
Wing root chord (CROOTW) 
Exposed wing area, 2 panels (SW) 


Wing mean eet chord (XMACW) 


Distance from nose 
Wing leading Sacae (XWING ) 


Wing thickness-to-chord ratio (TOVCH) 


Tail tip-to-chord ratio (XLAMT) 
Tail leading edge sweep (CLAMT) 
Tail span including body (BT) 
Tail root chord (CROOTT) 

Exposed tail area, 2 panels (ST) 


Tail mean geometric chord (XMACT) 


Distance from nose 
tail leading eee (XTAIL) 


Tail thickness-to-chord ratio (TOVCT) 


Altitude (HT) 

Body diameter (D) 

Missile length (XL) 

Nose length (XLNOSE) 

Distance from nose to CG (XCG) 
Reference area (AREA) 
Reference length (XREF) 

Engine code (ENGINE) 


Inlet code (ENLET) 


Boat tail angle (BETA) 

Base diameter (DBASE) 
Nozzle exit diameter (DJET) 
Boat tail length (XLABOD) 
Protuberance drag (CDPROT) 


If comparing with experimental 
values, Reynolds number (REFT) 


ot 


Feet 
Feet 
Square feet 
Feet 


Feet 


Degrees 
Feet 

Feet 

Square feet 
Feet 


Feet 


Feet 
Feet 
Feet 
Feet 
Feet 
Square feet 
Peet 


0.0-Turbofan 
1.0-Rocket 


Q0.0-Flush 
1.0-Extended 


Degrees 

Feet 

reet 

Feet 

(Coefficient value) 


(Dimensionless) 





12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problen. 

13. To raceive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 


press [ENTER] . 
i; Switch. 


15. Turn the terminal off with the red 


C. EXAMPLE PROBLEMS 
Table (V-5) identifies the output variables as they 
appear in the output tables. 
1. Example V-A. Tail control missile 
Figure (V¥-1) illustrates the missile used in this 
example. The dimensions for this missile and cther input 
parameters are contained in Table (¥-1). The output is 
Shown in Table (V-2) and Pigure (V-2). 
2.e Example V-It. Canard control missile 
Figure (V¥-3) illustrates the conard configuration 
Missile used in this example. The input data is contained 
in Tatle (V-3). The output is displayed in Table (¥-4) and 
Figure (V-4). 


ie 
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Pigure (¥-1). Tail control missile as used in Example V-A 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


There are many topics which may be the subjects of fol- 
low on work contained within this thesis. Although the four 
programs have keen installed on the IBM 370 computer system, 
these four alone do not fully satisfy the original goal of 
this work: Provide a computer supplement to the Tactical 
Missile Conceptual Design handbook. Numerous additional 
focal algorithms are utilized in the design handbook which 
deserve the attention of a programmer. Of immediate inter- 
est are the areas concerning radar or infrared guidance 
systems, baseline configuration modelling and weight distri- 
bution, and initial control and lifting surface design. 

Each of these topics can be programmed to provide missile 
design students interactive learning tools when coupled with 
the design handbook. 

The most urgent follow on work to this thesis is the 
restoration of the program LAERO1 to a reliable, useful pro- 
gram. The program was modified and set up on the IBM 370 
cOmputer system during the period immediately following the 
System's installation at the Naval Postgraduate School. As 
could be expected, the computer suffered many and varied 
grtcwing pains in its early life. As a result of this, or of 
the human manipulation expanding the capability of the pro- 
gram, the effectiveness of LAERO1 was substantially reduced. 

Work involving the other three programs would involve 
Slaply expanding their capabilities. The trajectory models 
program, LPATH, presently considers only two guidance laws: 
line-of-sight guidance and proportional navigation guidance. 
Other guidance laws which can easily be included in the pro- 
gtam include pursuit, beam rider, and combinations of diffe- 
rent laws. It might also prove useful to be able to 
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simulate the entire missile trajectory but still only output 
the terminal phase of the encounter. Another option would be 
tc provide the target with a controlled trajectory instead 
cf the constant acceleration condition now imposed. 

The propulsion sizing program, LPROP, should have the 
various nozzle options incorporated into the program so that 
it isn't necessary to manually juggle the program output. 
Other booster-Sustainer grain configurations could be 
explored, such as the booster grain being cast within the 
ccre of the sustainer, or even a motor with only a single 
grain. Another suggestion for the convenience of the proe- 
gram users is to institute a shopping list of availlable 
propellants and their characteristics into the progran. 

A Similac list of available explosives and case materi- 
als could be put into the warhead sizing program, LBOMB. 
These shopping lists would provide ready access to reference 
information and, at the same time, decrease the number of 
data values to be manually input into the computer. Since 
the current program is limited to cylindrical warheads, an 
area of expansion would be the flexibility of warhead 
styles, such as curved, shaped charge, continuous rod, etc. 

At present, the programs are somewhat hindered by the 
mechanics involved in producing the printouts and the plots. 
Due to the results of tailoring a program to be interactive, 
often it is required to completely exit a program before 
Gutput can ke received. Subsequent design iterations 
Leguire re-entry into the program, which produces a certain 
justifyable annoyance to the user. Additionally, the compu- 
ter center has instituted a new policy of cancelling any 
jobs with duplicate job names, which can be severely irri- 
tating and cumbersome to the persons running the plot rout- 
ines contained within this thesis. However, the computer 
center has developed procedures which have the potential to 
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alleviate both of these problems. According to Volume 13, 
Number 4 of the Computer Center Newsletter, CMS commands can 
now be invoked from within a FORTRAN program. fThe print and 
plot operations presently contained within executive rout- 
ines can now be placed directly within the source programs. 
This will remove all current restrictions placed on the nun- 
bers of printouts received per sesSion and will label each 
plot with the user's job name and not the project's job 


haze. 
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APPENDIX A. IRAJECTORY MODELS PROGRAM FLOWCHARY?Y 


MAIN PROGRAM 


NRUN=0 
[SAVE=1 






Input data: TPITLE,DEL,NAV,LOS 
IRT, TSP, TAC, TAL ,MSP,IRM 





[IPRO=1 
[TLOS=1 


Playback input data to terminal 
FOr Vent ication 


TTYI=LOS+TAL 
MTH=LOS+M AL 


OUVO0U0 
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ILOS=0 
[PRO=1 


ee. 


ee ean ge RT) 
MAL=A2 
MTH=MAL+LOS 


T2=TH 
sane 
[354 | 


=0 


KM (I) =0 
_ YM rf 
CT (I) =(a 1-82) =San(P1) 





Print 2rror message 
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“777 a5254c 7438 
P3=TSP*sin(A1) /R1 


D6=R4*TSP*sin (AT) 
D7=R2*TSP*A 3*coOsS (Al) 
Taree acne aD} 
DI=R1I*MSP*cos (A2 


XM (I) =XM (I-1) *+DEL*MSP*cos (T2 
YM (I) =YM (1-1) +DEL*MSP*sin (T2 
XT (I) =XT (I-1) *DEL*TSP*cos (T1 
YT =YT (I-1}) +DEL*TSP*sin (TI 











i -- 


[ Call OOTPUT 
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SP*cos 


M 

Mf 

iu ® 
TSP*sin 


+r ¥ 
H % -*—1 


~sCDB=NAVAMSPATAC cos (81) /8 (T-1) a 
D9=A (I-1) * (2#R2+NAV ¥*MSP¥COS (AZ)) /R (I-1) 


D2=(TSP*sin (A1) -MSP¥sin (A2)) /R (I-1) 
D6=TAC/TSP 


D7=NAV*D2 
P2=D6-D2 


on 





R (I) =R (I-1) +DEL*R2 
A{I) =A (I-1) +DEL*A3 
P1=P1+DEL*D2 
Ti=f1+DEL*D6 
T2=[2+DEL*D7 
A1=A1+DEL*P2 
AZ=A2+DEL*P3 









Do you want to run the 
Same problem again? 


Do you want +o run \_™ 
a new problen? 
No 
[ wRon=99 | 


re 
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APPENDIX BB. TRAJECTORY MODELS PROGRAM LISTING 


Following the next page is the program listing for the 
Trajectory Models program. It consists of three main divi- 
sions; the executiv2 routines, the FORTRAN IV computational 
program, and the PORTRAN IV plotting program. The #xecutive 
routines are used to establish required file definitions, 
initiate operation of the computational program, supervise 
the transfer of data to the plotting program, and provide 
Operational information to the program user at appropriate 
times. 

The computational program, LPATH FORTRAN, consists of 
four subprogram divisions. The MAIN program acc2pts the 
input data and performs th2 calculations for the line of 
Sight and the proportional navigation problems. Subroutine 
SCREEN is used to prompt the user to clear the terminal 
screen for proper positioning of the output. Subroutine 
SORT determines the largest missile acceleration value and 
the value ranges of the X and Y position coordinates for 
piotting reference. Finally, subroutine OUTPUT formats the 
Calculated data and provides it to tne user, the printer 
file, and the plot program data file. 

The plot progran, PATHPLOI FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
Minal. No cards ne2d to b= punched or read. Attention must 
be given to the correct JCL accounting data in the first 
line of the program; thos? shown in the listing ar2 for 
illustrative purposes only. fhis program receives the data 
from the computational program and produces a singls geo- 
graphic chart of the encounter in the encounter plane. Mul- 
tiple problems, up to nine, will overlay on the single 
emart. 
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fmeues LPATH EXEC A NAVAL POSTGRADUATES SCHOOL 
FILEDEF 08 DISK LPATH OUTPUT AO 

FILEDEF 07 DISK PATHPLOT DATA AQ 

&BEGTYPE 

YOU WILL HAVE THE OPTION [O OBTAIN BOTH A HARDCOPY PRINTODUT 
meoeA VERSATEC PLOT OF UP TO NINE ENCOUNTER GEOMETRIES. 

THE PLOT IS A SINGLE FRAME WITH ALL NINE GEOMETRIES 
SUPERIMPOSED ON ONE ANOTHER. THE HARDCOPY PRINTOUT IS 
IDENTICAL IN FORMAT TO THE TERMINAL OUTPUT. 


& END 

LOAD LPATH 

START 

&SBEGTYPE 

meee. oN A HARDCOPY PRINTOUT OF THE RESULTS, TYPE AND 


lpathpr 
THE OUTPUT WILL BE PRINTED ON THE VM PRINTER IN ROOM 140 
meeewliLL BE IDENTIFIED B33 YOUR USER NUMBER. EL UsuALLeyY 
MeouLRES SEVERAL MINUTES [TO OBTAIN THE PRINTOOT. 
TO OBTAIN THE PLOT OF YOUR ENCOUNTERS, TYPE AND ENTER: 
lpathpl 
meee LOT WILL BE PRINTED IN THE COMPUTER ROOM AND WILL BE 
PLACED ON TOP OF THE ALPHABETIZED OUTPUT FILE IN ROOM 140. 
feeerul BE LI DENTIFLTED BY fHE JCL JOB NAME “PATHPLIT" AND 
USUALLY REQUIRES MANY MINUTES Le HOURS!) TO BE 
MeeeeNeED. NOTE...-IF MANY USEKS ARE REQUESTING PLOTS 
PLL ANEOUSLY, THE COMPUTER CENTER PERSONNEL WILL CANCEL 
"EXCESS" JOBS OSINS Bie “AWG DDENTEPLER. 
& END 
Meoes LPATHPL EXEC A NAVAL POSTGRADUATE SCHOOL 
COPY LPATH PLOT A PATHPLOT DATA A PLOT FORTRAN A 
EXEC SUBMIT PLOT FORTRAN A 
ERASE PLOT FORTRAN A 
Moe, LPATHPR EXEC A NAVAL POSTSRADUATE SCHOOL 


PRINT LPATH OUTPUT (LINECOUN 70 
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APPENDIX C. WARHEAD DESIGN PROGRAM FLOWCHART 
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"RHO=.09237 698 (exp (- - 0000 3254ALT) = .0172 
*exp (.0 eee ee woe) 


NO 
ALIT< 36000? 


YES 


NO 
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epecs= (CORRS SPI /ERALB(T))* 2* RW*RW*CTH-RW*CTH*CTH) 
PRLEN= (FRMLB(1) / (CTH*CDEN) ) **.5 
HITS=FRAGS#A 

Return ] 


ES 





APPENDIX D. WARHEAD DESIGN PROGRAM LISTING 


This program has two major sections; the executive rout- 
ine and the FORTRAN IV computational program. The executive 
routine establishes the required file definitions and initi- 
ates operation of the computational progran. 

The computational program, LBOMB FORTRAN, consists of 
four subprogram divisions. The MAIN program accepts the 
input data, calculates the atmospheric characteristics, and 
formats and displays the output to the user and sends it to 
the printer file. Subroutine VEL calculates the initial 
velocity required to propel a given mass a specified dis- 
tance through the atmosphere with a particular residual vel- 
ocity remaining. It also determines the ballistic limit 
velocities for the situation. Subroutine SIZER siz2s the 
warhead for a given Pd value. It also produces the charge- 
to-mass ratio, the number of fragments, the fragnent size, 
and the average number of hits received by the target. 
Subroutine SCREEN prompts the user to clear the terminal 


screen for proper positioning of the displayed data. 
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SEND 


SEND 
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APPENDIX E. PROPULSION SIZING PROGRAM FLOWCHART 


| MAIN PROGRAM | 


Input IMO 
DRAGS,&, 








tr3jr3 









NO 


Input ALTBN,ER,DB,YIELD,DENSC 


Input ALTBN 
ALTSN, 





NCSC 





DV=VBP-VB 
Tats 
DRAGB= ((VBP+VBTI) * 


[ WPB=TBB*DM 


eee eee 
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TIB=WPB*ISPB 
FB=TIB/TBB 
B= (DV/ 6*788)) ¢sin (LB) 
WN=WL* (1-exp ( -WPB/(FB-DRAGB) ) *B) ) 












PCB=PAB 
PEB=PAB 
AR=1 


Call PCRAT 


ae eae 
Call ARAT 


— ent? = ap ap ef-s <= => =p 








Is 
IMOTO R=0? 
or 
DEB> (0.95*DB) ? 
or 
PCB<1000? 





ND 


PLIM=PLIM*.99 
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2*LB 
Call PRATIO 
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BRB 
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Call PRATIO 


DRAGS*1.1+ (WL-TWB) *sin (ELS) 


ot 
L ESSDERGST IA 1¢ (0 cPED) tata (ELS) 








PCS=PAS 
PES=PAS 
AR=1 





Call ARAT 






DES=0. 95*D? 






PLEM=PLIN=.99 
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Set all sustainer values 
C5 2Zero. 
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L=LB+tLS +NLB+tNLS 
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Call OUTPUT 
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( 
S 
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@ oo - 
NO 


PC=PT* 1.01 






(PE/PA) >0.4? 1 


(PE/PA) <2.5? 


Return | 
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APPENDIX F. PROPULSION SIZINS PROGRAM LISTING 


Following this page is the program listing for the Pro- 
pulsion Sizing Program. It has two segments; the 2xecutive 
routine and the FORTRAN IV computational program. The exe- 
cutive routine establishes the required file definitions and 
initiates operation of the computational program. 

The computational program, LPROP PORTRAN, consists of 
Six subprogram divisions. The MAIN program accepts the 
input data from the user and performs the guiding calcula- 
tions for the booster and sustainer motors, Subroutine 
PRATIO determines the ambient pressures at the design alti- 
tudes. Subroutine PCRAT defines the optimum chamber pres- 
Sure to ambient pressure ratid with respect to the case 
Material properties. Subroutine ARAT solves for the area 
ratio of the nozzle and tries to size the nozzle to the mis- 
Sile diameter by varying the chamber pressure, characteris- 
tic velocity, and thrust coefficient. Subroutine SCREEN is 
used to prompt the user to clear the terminal screen for 
proper positioning of the sutput. And subroutine OUTPUT 
formats the calculated solutions and provides them to the 


user and to the printer file, if so directed by the user. 


ee 





NAVAL POSTGRADUATE SCHOOL 


A 


: DEF 08 DISK LPROP OU ROUT AQ 


ELLE 
&BEGTYPE 


EXEC 


LPROP 


FILE 


(RECFM VA BLOCK 133 PERM 


LOAD LPROP 


START 
& BEGTYPE 


& END 


& END 
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APPENDIX G. AERIJDYNAMIC COEFFICIENTS PROGRAM LISTING 


This program is divided into three major subdivisions; 
the executive routines, thea FORTRAN II/IV computational pro- 
gram, and the FORTRAN IV plotting program. The executive 
routines are used to establish the required file defini- 
tions, initiate »speration of the computational program, 
supervise the transfer of data to the plotting program, and 
provide operational information to the user at appropriate 
times. 

The computational program, LAERO! PORTRAN, consists of 
five subprogram divisions. The MAIN program accepts the 
input data, conducts calc lations not done by other subrout- 
ines, formats the calculated data and provides the sutput 
data to the user, the printer file, and tha plot program 
data file. Subroutine GEOSUB calculates the missile wetted 
area and the Reynolds number per foot based on the flight 
altitude. Subroutine CLASJB calculates the aerodynamic sur- 
Face lift-curve slopes. Subroutine CATSUB calculatas center 
of pressure locations, cross-flow drag coefficients, and 
interference factors. Subroutine SCREEN prompts the user to 
clear the terminal screen for proper positioning of the out- 
put. 

The plot progran, AEROPLOT FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
Minal. No cards need to ba panched or read. However, 
attention must be given to the correct JCL accounting data 
in the first line of the program; those shown in the list- 
ing are for illustrative purposes only. This progran 
receives the data from the computational program and pro- 
duces a group of six charts for @ach Mach number entered. 
The plots represent the relationships of Cl to alpha, Cn, 
Eiyeca, Ca, and Cdi. The program can produce up t> 24 sets 
Of plots for a single run. 


£53 





PELE: LAEROT EXEC A NAVAL POSTGRADUATE SCHOOL 


FILEDEF 08 DISK LAERO1 OUPOr AD (RECFM VA BLOCK 133 PERM 
oe e. DISK LAERO1 PLOT AO (RECFM VA BLOCK 80 PERM 
YOU WILL HAVE THE OPPORTUNITY TO OBIAIN BOTH A HARDCOPY 
PRINTOUT AND A SET OF PLOTS FOR ONE SET OF INPUT PARAMETERS 
EACH TIME YOU ENTER THIS PROGRAM. IHE PROGRAM MAY BE 
RE-RON CONTINUOUSLY AND YOU WILL HAVE THE OPTION TO CHANGE 
INPUT PARAMETERS FOR EACH SUCCESSIVE RUN, BUT YOU CAN 
OBTAIN THE PRINTOOT AND PLOTS PERTAINING [TO THE LAST RUN 
ONLY. IF ADDITIONAL OUTPIJT IS REQUIRED, RE-ENTER THE 
PROGRAM. 
& END 
LOAD LAERO1 
START 
& BEGTYPE 
TO OBTAIN THE HARDCOPY PRINTOUT OF THE DATA TABLES, TYPE 
AND ENTER: 
LAEROIPR 
poeple THE VERSATEC PLOT OF THE TABULAR DATA, IYPE AND 
LAEROIPL 
& END 
FILE: LAEROWPR EXEC A NAVAL POSTGRADUATE SCHOOL 
PRINT LAERO1 OUTPOT (LINEZIIUN 70 
&SBEGTYPE 
PaeeOUTPUT WILL BE PRINTEO ON THE VMS PRINTER IN ROOM 1740 
AND WILL BE IDENTIFIED WITH YOUR USER NUMBER AND LAST NAME. 
> acer! Pesos v ona Lea ENU TES To OBTAIN THE PRINTOUT. 
FILE: LAEROUPL EXEC A NAVAL POSTGRADUATE SCHOIOL 
COPY LAERO1 PLOT A LAERO PLOTDATA A PLOT LAERO! A 
EXEC SUBMIT PLOT LAERO1 A 
ERASE PLOT LAERO1 A 
&SBEGTYPE 
THE PLOT WILL BE DRAWN IN THE COMPUTER ROOM AND PLACED OVER 
THE ALPHABETIZED OUTPUT FILE IN ROOM 140. IT WILL BE 
PBEN TIFIED BY ROMO AND USUALLY REQUIRES MANY MINUTES 
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